Study of the reaction pp (pep from 1.1 to 2.0 GeV/c 
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A study has been performed of the reaction pp 4K^ using in-fiight antiprotons from 1.1 to 2.0 
GeV/c incident momentum interacting with a hydrogen jet target. The reaction is dominated by 
the production of a pair of (j> mesons. The pp (jxj) cross section rises sharply above threshold and 
then falls continuously as a function of increasing antiproton momentum. The overall magnitude of 
the cross section exceeds expectations from a simple application of the OZI rule by two orders of 
magnitude. In a fine scan around the 5//j(2230) resonance, no structure is observed. A limit is set 
for the double branching ratio B(^ — > pp) x B(^ <j)(j)) < 6 x 10^^ for a spin 2 resonance of M = 
2.235 GeV and T = 15 MeV. 



I. INTRODUCTION 

Beyond the common famihes of baryonic {qqq) and 
mesonic {qq) states, glueballs {gg or ggg), hybrids {qqg) 
and multi-quark systems {e.g., qqqq) are predicted to ex- 
ist based on the current understanding of QCD and 
QCD-inspired phenomenological models of hadrons. The 
question of the existence of these unusual hadronic bound 
systems is fundamental and should lead to a broader 
understanding of the behavior of QCD in the non- 
perturbative region. The experiment described below in- 
vestigates the existence of gluonic hadrons in the mass 
range between 2.1 and 2.4 GeV/c^. The technique re- 
quires a state, if found, to couple appreciably to both 
antiproton-proton (the entrance channel) and 4>4' (the 
exit channel) . The parameter space explored is of partic- 



ular interest in a search for the tensor glueball. 

Gluonic hadrons are expected to populate the low mass 
region of the hadron spectrum together with normal qq 
states. Recent lattice QCD calculations locate the scalar 
(J-^*-^ = 0"'""'") glueball mass in the range of 1.5 to 1.7 
GeV [§. The tensor [J^^ = 2++) is the next low-lying 
glueball and is found to be approximately 1.5 times more 
massive than the scalar, thus it is postulated that it might 
be in the region explored by this experiment. 

Several experimental approaches have been used to 
search for glueballs. The common theme is to choose a 
process where a kinematical or dynamical filter strongly 
reduces the large production of standard qq mesons. Re- 
actions intensively explored are peripheral hadron inter- 
actions 1^ , J/i/) radiative decays [|j , central production 
Q , 77 collisions ||] and pp annihilations Q| . 
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The study of the pp (j)^ reaction is motivated by 
the expectation that, with essentiahy no common va- 
lence quarks between the initial and final states, the re- 
action is forbidden in first order. Indeed, a strict appUca- 
tion of the empirical Okubo-Zweig-Iizuka (OZI) rule |^] 
suggests a cross section for pp — > (fxj) of approximately 
10 nb. In a previous publication we reported the 
cross section at 1.4 GeV/c incident antiproton momen- 
tum to be 3.7 ± 0.8 /Lib and thus two orders of magnitude 
greater than estimated. The large pp — > cross sec- 
tion has been interpreted as coming from two-step pro- 
cesses [p|JIO|], from intermediate glue and/or from 
the intrinsic strangeness in the nucleons |12(|. 

Here we review the appropriateness of the energy range 
of our study with respect to the possibility of exciting a 
gluonic resonance. As evidence, consider the current pic- 
ture of candidate glueball states. The /o(1500), discov- 
ered by the Crystal Barrel experiment at the CERN Low- 
Energy Antiproton Ring (LEAR), has a width of F « 120 
MeV and has been observed to decay to tt^tt", 7777 [ p^ , 
7777' jlj] and KK [|l^ following pp annihilation at rest. 
The existence of this state has been confirmed in in-flight 
pp annihilations by experiment E760 at Fermilab JTot . It 
has been interpreted as the leading scalar glueball candi- 
date. Alternate explanations feature strong gluonic con- 
tent; the state could be a mixture of a glueball and nearby 
ordinary qq systems [|l^. If identified with the scalar 
glueball, the implication from the lattice calculations is 
that the tensor mass should be in the 2.1 to 2.4 GeV 
range where candidate glueball states have already been 
identified. These states include the ^//j(2230) which 
is observed in radiative J/ip decays to KK |l^ and 
also in decays to tttt and pp |19|. The width of this 
state (F « 30 ± 20 MeV) is unexpectedly small, and has 
stimulated considerable interest. The spin-parity is un- 
determined but is limited to (even) + . This channel is 
of interest to the present experiment since it couples to 
pp. In another study |20[| from this experiment, strin- 
gent limits were established on the double branching ra- 
tio B(^ pp) X B(^ KgKg) for the case of a narrow 
(i.e., F < 30 MeV) state. 

Resonances in the (jxj) system have been observed in 
J /ip 700 by Mark III and DM2. Each group saw a 
wide bump above (jycf) threshold peaking near 2.2 GeV and 
having a width of approximately 100 MeV. The preferred 
assignment is 0^, however 2+ contributions cannot be 
excluded |2|. 

The (jxj) system has been studied in exclusive -K^p 
p2| , p3t and inclusive it~ Be induced interactions pi] . 
This reaction is expected to be highly suppressed based 
on the above OZI argument related to disconnected 
quark-line diagrams. However the suppression is not ob- 
served near threshold, a fact which could be explained by 
the production of intermediate gluonic resonances [ pT| . 
A partial wave analysis (PWA) of the exclusive data re- 
vealed the presence of three interfering J^'-^ = 2++ res- 



onances, the /2(2010),/2(2300) and /2(2340), originally 
called gx states. The inclusive data show evidence for 
two structures with parameters compatible with those 
measured for the qt states but a PWA was not possible. 

We report below on the evaluation of 58 independent 
measurements of the reaction pp — > AK^ from which we 
determine the cross sections for pp — > 00, pp — > <j>K'^ K~ 
and pp — » AK^. A threshold enhancement is evident in 
the 00 cross section data. In the narrow region around 
the ^//j(2230) resonance, 17 measurements of the 00 
cross section are reported in narrow energy steps covering 
a range of 45 MeV. Limits are set on the non-observation 
of this state. 



II. THE EXPERIMENT 

The JETSET (PS202) || experiment is designed to 
measure the reaction 



pp 



K+K-K+K 



(1) 



as a function of incident antiproton momentum from 1 . 1 
GeV/c to the highest momenta available at LEAR (2.0 
GeV/c). The experiment makes use of the LEAR an- 
tiproton beam incident on an internal target provided by 
a hydrogen-cluster jet system. To obtain a measure of the 
absolute (and relative) luminosity at different antiproton 
momenta, pp elastic scattering events were recorded in 
parallel using a special dedicated trigger. 

One key feature of reaction (1) in the energy range 
covered by this experiment is the fact that the outgoing 
kaons are constrained to a forward cone below 60°. Dom- 
inant in pp annihilation are the production of charged 
and neutral pious. These unwanted background reac- 
tions produce relatively fast charged particles and pho- 
tons spread over a much larger angular range. The design 
of the detector and electronic trigger takes advantage of 
these facts by providing good charged-particle tracking 
and particle identification (PID) for forward angles, fast 
charged-particle multiplicity for triggering, and a large- 
acceptance photon rejection system. The philosophy fol- 
lowed to record and identify 4,K^ events is to trigger only 
on the proper multiplicity pattern. In the offline analysis, 
events are required to have four well-determined tracks 
in the right kinematical range. Finally, these tracks are 
associated with hits in the PID detectors to determine 
compatibility with identiflcation as kaons from the AK 
final state. 

A schematic view of the apparatus is shown in Fig. 
Its basic structure is a hydrogen-cluster jet target (not 
shown) with a density p = 4 x 10^^ atoms/cm^, sur- 
rounded by a compact detector. The jet is oriented hori- 
zontally and intersects the circulating antiproton beam at 
right angles. Around the interaction volume the LEAR 
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ring is equipped with an oval vacuum chamber (0.03 ra- 
diation lengths thick) having horizontal and vertical half- 
axes of 78 and 38 mm, respectively. This limits the 
geometrical acceptance of the detector to polar angles 
0iab > 7°, with complete coverage of azimuthal angles 
only for 9iab > 15°. 

The overall structure of the detector includes a "for- 
ward" and a "barrel" sector with polar angular coverage 
from 7° —45° and 45° — 135°, respectively. Immediately 
outside of the beam pipe is an "inner" scintillation ho- 
doscope array segmented in azimuth into 40 (forward) 
and 20 (barrel) components. These are used in conjunc- 
tion with a triple-layer "outer" hodoscope (144 total el- 
ements) |27| to form multiplicity patterns for use in the 
trigger. In the reaction pp AK^ at least three charged 
tracks are required within a forward 45° cone around the 
beam axis. The fourth kaon has a maximum polar angle 
of 60°. The inner scintillator array extends to this angle, 
while the barrel outer array includes angles up to 135° 
to veto events with large angle tracks. 

Between these scintillator arrays are cylindrical drift 
chambers ("straws") which are used for tracking [ p8| . 
They are arranged in 12 planes with alternating hori- 
zontal and vertical orientation in the forward region. In 
the barrel, 1400 straws are formed into a tightly packed 
bundle aligned with and surrounding the beam pipe. 

Three separate devices are incorporated with the pur- 
pose of providing PID information. These include 3500 
silicon-pad counters which measure energy loss. They 
are particularly effective at the low end of the expected 
(3 range. A hodoscope consisting of 24 forward and 24 
barrel threshold Cerenkov counters is positioned just in- 
side of the outer scintillator array. These counters are 
used to limit the number of "fast" particles that are ac- 
cepted at the trigger level (typically < 2) and, for those 
events which induce light in the Cerenkov counters, as 
a measure of the /3 of that passing particle. Finally, a 
RICH detector, having a quartz radiator {(3th — 0.64), 
was introduced into the forward detector compartment 
following the first year of data taking It is capable 
of providing a precise measurement of the (3 for typically 
one of the forward-going tracks in each event. Outside 
of all of these components an electromagnetic calorime- 
ter made from lead and scintillating fibers (Pb/SciFi) 
is placed. In the forward region the shower counter is 
divided into 300 towers. Each tower is 12.5 radiation 
lengths thick and has individual phototube readout . 
In the barrel, 24 trapezoidal bars of Pb/SciFi form a 
hermetic cylinder with six radiation lengths thickness at 
normal incidence. The barrel calorimeter, in conjunc- 
tion with the outer barrel scintillator array, has a special 
electronic pattern unit enabling the fast identification of 
isolated high-energy photons. The forward calorimeter is 
passive, used only in the offline analysis for evidence of 
high-energy gammas in the event. 

The online trigger was based on a charged-particle mul- 




RICH Inner Scintillators 

FIG. 1. Layout of the JETSET/PS202 experiment, show- 



ing the major detector components. A typical event with 
three forward and one barrel track is shown. 

tiplicity of four (with no more than one track in the bar- 
rel region) coupled with a crude transverse momentum 
balance requirement based on the forward straws mul- 
tiplicity. Events with more than two of the threshold 
Cerenkov counters responding or with high-energy gam- 
mas in the barrel calorimeter were vetoed. These trigger 
conditions were tested using special runs in which such 
conditions were relaxed in order to determine the sen- 
sitivity of the detector acceptance to the trigger. For 
events meeting the trigger conditions, information was 
recorded from multiple parallel VME-based subsystems 
onto a common data stream pl| which was recorded on 
magnetic tape. Approximately 500 million events were 
recorded in eight running periods spanning the time from 
April 1991 through September 1994. 

III. DATA SELECTION 

A. Event reconstruction 

The offline analysis proceeds as follows. Since the de- 
tector is non-magnetic, only the directions of the four 
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FIG. 2. Ai5 (see text) distributions for (a) Monte Carlo simulation of the reaction pp AK^ at 1.5 GeV/c; (b) Raw 
four-prong data for different momenta ranges: Light line: 1.2 GeV/c, grey line: 1.65 GeV/c, black line: 2.0 GeV/c; (c) 
Application of the dE/dx criterion only to data in the 1.4-1.45 GeV/c region; (d) Adding the threshold Cerenkov compatibility 
condition to this data; (e) Adding the criterion from the RICH and from the 7 veto; (f) Final AS distribution for all the data 
taken in the experiment. The shaded area shows the cut used to select the reaction. 



outgoing particles in the reaction 

pp^ K+R'-K+R'- (2) 

are measured. The four momenta therefore are obtained 
by solving the four equations of energy-momentum con- 
servation. 

The conservation of three-momentum from the well- 
defined initial state (e.g., Ap/p « lO^'^ for the antipro- 
ton) allows the magnitudes of the momenta for the four 
kaons in the final state to be expressed as linear functions 
of one unknown parameter, which can be taken to be the 
momentum, p^, of the fourth kaon. The total energy in 
the laboratory of the final state E = J2 \/Pi ~^ ™/f 
then be written as a function of this one parameter p^. 
The function E{p4) has a roughly hyperbolic shape that 
diverges as J34 ±00 and possesses a single minimum. 
We define /S.E — E{p4) — Eq as the minimum of the func- 
tion E{p4) minus the total energy Eo of the initial state in 
the laboratory system. This variable AE is constrained 
to AE < for events from reaction (1), with a singu- 
larity appearing in the available phase-space at AE — 0. 
The width of this spike is determined mainly by the res- 
olution of the straw chamber tracking system which is 
well studied. The resulting resolution on the momentum 
determination is 8% at a typical kaon momentum of 0.5 
GeV/c. 



A Monte Carlo simulation of pure pp AK events, 
with the experimental resolution and all other material 
effects included, results in a AE distribution as shown in 
Fig. |^(a) . The distribution is asymmetric, and its width 
increases with increasing incident p momentum. 

The four-prong data, prior to application of any PID 
information, show a broad AE distribution with no obvi- 
ous structure. This is consistent with expectations since 
these events are dominated by background. The shape 
variation of such raw AE distributions over the span of 
incident momenta explored by our experiment is shown 
in Fig. |(b). 

B. Particle Identification 

In order to isolate the pp AK^ events from the 
very large background, use of different selection crite- 
ria based on the particle identification devices is made. 
First we require an energy loss (dE/dx) compatibility as 
measured in the silicon detectors. This is evaluated by 
forming a confidence level for each measurement, using a 
parametrization of the Landau distribution and integrat- 
ing the tail beyond the measured value. The confidence 
levels for each measurement are then combined into a 
dE/dx confidence level. The detectors were calibrated 
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FIG. 3. (a) Pulse height in ADC units for the barrel silicon 
pads vs. P for a sample of elastic pp events, (b) Pulse height 
in ADC units for the forward silicon pads vs. (3 from a sample 
of p or p from the reaction pp 
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FIG. 4. (a) Measured versus expected /3 for elastic events 
in the forward Cerenkov counters with the water radiator, 
(b) The fraction of tracks giving light in the Cerenkovs as 
a function of (3 for the freon radiator. Events of the type 
pp pp-K^-K~ are used. Light observed below threshold re- 
sults from the plexiglass walls of the counters. 



using a sample of well-defined tracks with known mo- 
menta. They come from fully identified pp elastic scat- 
tering events and from pp — > ppiT^TT^ events. This latter 
reaction was isolated almost background free by making 
use of the Cerenkov identification and by the requirement 
of a large energy deposit from the outgoing p in the for- 
ward electromagnetic calorimeter . A plot of the pulse 
height in the silicon pads vs. (3 is shown in Fig. ^(a) for 
the barrel region and in Fig. ^(b) for the forward region. 
The first plot makes use of elastic events while the second 
distribution has been obtained from pp —> pp7r+7r~ data. 
Application of the dE/dx compatibility requirement on 
a sample of the 1.4 to 1.45 GeV/c data results in a Ai? 
distribution as shown in Fig. |^(c). At this stage the 
pp — > 4:K^ events are already in evidence. 

The Cerenkov effect produces a response R{I3) in pho- 
toelectrons that is monotonically increasing for /3 > 
(3th- Two different radiators were used: liquid freon 
(FC72) |3^], in the very first data taking only, having 
a threshold at /3=0.79, and subsequently water having 
a threshold at 0.75. The result of the (3 measurement 
from the Cerenkov response events of the water radiator 
is compared in Fig. ||(a) to the expected [3 from elastic 
events. 

For the freon data, we show in Fig. ^(b) the "effi- 
ciency" versus (3 where we have defined efficiency as the 
ratio between the (3 distribution of the tracks crossing the 
Cerenkov counters which give light to the /3 distribution 
of all the tracks. This study is based on events of the 
type pp ppTT~^TT~ . A smooth increase as a function of 
(3 is seen which indicates the expected threshold behav- 
ior. The non-zero efficiency below threshold is caused by 
light produced in the plexiglass containment wall of the 
counters. 

The response function of the liquid Cerenkov counters 
is used to calculate the expected signal from a passing 
track of a given hypothesized momentum. The expected 



and measured responses are compared on the basis of 
Poisson statistics. The resulting confidence levels for the 
four tracks are then combined with the one obtained from 
the dE/dx information to form an overall particle iden- 
tification confidence level which is required to be at least 
5%. The effect of adding the Cerenkov information can 
be seen in Fig. ^(d). 

The response function of the RICH is determined by 
measurements of elastic events in which the forward- 
scattered antiproton of known momentum penetrates the 
detector. Details of this detector performance are found 
elsewhere 29 1 , however a representative example of the (3 
resolution is shown in Fig. ^. The difference between the 
expected and measured (3 distribution exhibits a Gaus- 
sian behavior with a resolution of /S.j3/ (3 ~ 2%. In ana- 
lyzing the response of the detector, we find that a large 
fraction of background events produces a RICH (3 mea- 
surement which is unphysically greater than 1. The most 
effective cut eliminates this background, yet retains the 
good events, by simply requiring that the [3 measured by 
the RICH be less than 0.9. 

One background reaction which frequently satisfies the 
trigger conditions of the experiment is 



pp —> ppn TT 



(3) 



which is particularly important for p momenta above 1.6 
GeV/c. Reaction (3) is identified by the same selection 
criteria as those applied to find pp AK^ events. Events 
so identified are removed from the final AK^ sample. 

At the final stage of event selection, those events with 
isolated 7's in the barrel region are removed. This elim- 
inates most of the events with multiple 7's or tt'^'s. The 
final /S.E distribution after applying all the above selec- 
tion criteria is shown in Fig. ^e) for the 1.4 to 1.45 GeV/c 
sample and in Fig. |^(f) for all of the data taken in the 
experiment combined. A clear signal, peaked at a IS.E 
of zero, is seen which is evidence for events of the type 
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FIG. 5. Expected minus measured /3 from the RICH de- 
tector for elastic events. 

pp AK^. The final selection of events is made by re- 
quiring /S.E < 20 MeV. This cut selects approximately 
32 000 events which are retained for further study. 

IV. ACCEPTANCE AND LUMINOSITY 

The experiment accumulated data at 58 incident p mo- 
mentum settings with different luminosities over a period 
of four years. However in the following analysis, measure- 
ments at the same or nearly the same momentum have 
been grouped together in order to reduce statistical and 
systematic errors. The distribution of luminosities as a 
function of the momentum are summarized in Table I. 

The luminosity was monitored via the known cross sec- 
tion of elastic pp scattering. Two independent methods 
were employed to measure it. In the first method, the 
rate of coincidence between opposite pixels formed in the 
outer trigger scintillator hodoscope was used to count 
elastic events near 90° in the cm. system. The second 
method made use of small microstrip detectors placed 
near 90° in the barrel to detect the recoil particle from 
low-angle elastic scattering |Q . The data have been nor- 
malized as a function of time by using minimum bias 
scalers which were found to be very reliable and stable 
over the lifetime of the experiment. The relative lumi- 
nosity error is 2%. 

The acceptance for the elastic and AK events was com- 
puted by a GEANT ^ based Monte-Carlo simulation of 
the physical detector and trigger conditions. The Monte- 
Carlo events were subjected to the same reconstruction 
criteria and analysis cuts as the real data. We generated 
Monte-Carlo events for each of the 58 momentum set- 
tings, taking into account the changes in the layout of 
the apparatus and trigger conditions. For the physics 
analysis, three reactions were generated: pp (jxj), 



pp (j)K'^K~ and pp K~ K~^K~ . The simula- 

tion assumed phase space distributions for all three re- 
actions. The distributions of the total acceptance as a 
function of the incident p momentum are shown as fitted 
polynomials in Fig. ||(a). 

The differential acceptance function for pp (j)(f) is 
greatly affected by the kinematics of the reaction and 
the layout of the detector. The distribution of cosQcm is 
shown in Fig. ^(b) for all the data. Here, 9cm is defined 
as the centre-of-mass scattering angle of the outgoing <j>. 
A strong depletion of events occurs in the forward region, 
due to loss of tracks in the beam pipe. On the other 
hand the acceptance in the x angle is almost flat as it 
can be seen in Fig. 1(c). Here x is defined as the angle 
formed by the two 0's decay planes, in the 00 centre of 
mass system. Fig. ^(c) shows the distribution of the x 
angle for the data in the 1.26-1.65 GeV/c region (points 
with error bars), together with the distribution of x for 
Monte-Carlo phase space generated data folded with the 
acceptance of the apparatus. This missing forward angle 
acceptance in 9cm affects our ability to extract total 00 
cross sections as described later. 

Systematic errors on the cross section measurements 
have been estimated using the known cross section for 
pp — > pp'K'^'K~ [ ^2| whose determination was found to 
be in good agreement with previous measurements. We 
also used the fact that the cross sections at several en- 
ergy points were measured several times under different 
experimental conditions. We estimate the overall uncer- 
tainty on the absolute scale of the 00 cross section to be 
20 %. 

V. THE DATA 

For the selected pp event candidates 

we show in Fig. in Fig. ^(a) and in Fig. |8|(b) the scatter 
plot of the invariant masses m{K^K4) vs. m{KiK2)- 
Fig. [t] shows all the available data as a surface plot, where 
a strong accumulation of events can be observed at the 
nominal 00 position. Fig. ^(a) and (b) illustrate the same 
data under the form of a scatter diagram for two different 
regions of incident p momentum, i.e. below and above 1.5 
GeV/c respectively. Notice that, due to the absence of 
charge information, three combinations per event enter 
in these plots. We observe in these distributions a strong 
enhancement at the position of the 00. In the higher 
momentum region we also observe horizontal and vertical 
bands which indicate the presence of the pp — > (pK'^ K~ 
final state. Monte-Carlo simulations confirm that the 
diagonal bands represent the reflection of the 00 peak 
due to the multiple combinations. Selecting one 0, i.e. 
requiring one •m{Ky,K4) combination to lie in the region 
1.00-1.04 GeV/c^ and plotting the opposite combination 
m{KiK2), we obtain the distributions shown in Fig. 0(c) 
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FIG. 6. (a) Acceptance for the three reactions pp — > (j>cj>, pp 4>K'^K~ and pp AK^ . Distribution of (b) cosScm and 
(c) X for the reaction pp — > 00 in the momentum range between 1.26 and 1.65 GeV/c. The data are represented with error 
bars. In (c) the solid line represents the distribution expected for phase space. 




FIG. 7. m{Ks, K4) vs. m{Ki, K2) (three combinations 
per event) for all the events. 



and (d) where a clean peak is visible. The structure is 
well centered at the nominal d> mass. 



VI. CHANNEL LIKELIHOOD FIT 

In order to separate the (fxp cross section from the 
(j)K^K~ final state and from the and other back- 
ground mixture, the channel likelihood technique |Q was 
used. The method performs a maximum likelihood fit to 
the data using three amplitudes: 00, cpK^K^ and phase 
space (which is a mixture of AK and background and can- 



not be separated at this stage). The different channels 
have been described by the following amplitudes: 

3 

00 : A^^ = ^ Bi{mKK) x Bj{mKK), 



(j)KK : A^KK 



E 



BiirriKK) 



The likelihood function employed is the following: 



A. 



, A^KK /, N 

+ X^KK-^ 1" (1 - X^cf, - X^KKj- 

^rpKK 



In the above expressions a;^^ and x,pKK represent the 
fractions of the 00 and (j)K~^K~ channels respectively, 
Bi{mKK) represents the lineshape functions, I^^ and 
I<pKK represent the corresponding normalization inte- 
grals for the two amplitudes describing the two channels 
respectively. 

The width is a narrow resonance (r=4.4 MeV) there- 
fore the observed width in this experiment is dominated 
by the detector resolution which we found to be non- 
Gaussian. In order to describe the lineshape we made 
use of the Monte-Carlo simulations. The best representa- 
tion of this lineshape is obtained by means of a relativistic 
spin Breit-Wigner form whose full width varies from 11 
to 15 MeV as the incident p momentum increases from 
1.1 to 2.0 GeV/c. 

The normalization of the amplitudes was computed nu- 
merically using the Monte-Carlo simulations of the reac- 
tion pp AK. The resulting integrals were then fit to 
polynomials as functions of the incident p momentum. 
Changes in the layout of the apparatus and trigger con- 
ditions have little effect on the values of these integrals as 
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FIG. 8. m(ii'3,if4) vs. m{Ki, K2) (three combinations 
per event) for a) p incident momentum below 1.5 GeV/c and 
b) above 1.5 GeV/c. c), d) mass projections of a) and b) in 
the band, i.e. plots of m{KiK2) for 1.00 < m{K3, K4,) < 
1.04 GeV/c^ 



seen in Fig. ^ where these integrals and the fitted poly- 
nomials are shown. 

The channel likelihood fit gives the fractions of events 
as well as a probability for each event to belong to one of 
the three hypotheses. The results of the fits are visually 
displayed in Fig. Here we show the combinatorial 
m{KK) distribution (six entries per event) for incident 
p momentum below and above 1.5 GeV/c. The white re- 
gion represent the (fxj) contribution, the grey region shows 
the (j)K^K~ contribution while the black region shows 
the phase space (background -I- 4K) contribution. We 
notice that no peak has been left in the phase space 
distribution, indicative of a successful fit. We also notice 
the strong increase of the (j)K~^K~ contribution in the 
higher momentum region. The analysis gives a total of 
approximately 11 400 (jfcf) events. 

The channel likelihood method is able to separate the 
three (j)(f>, (f)K^K~ and (background -I- AK) contribu- 
tions. In order to separate the AK contribution from the 
background, the AE' distributions as discussed in section 
III. A. For this purpose, in order to reduce the errors, a 
further compression of the data was performed, grouping 
them in nine slices of incident p momentum. The /S.E dis- 
tributions for all intervals are shown in Fig. A clean 
peak at /S.E = over some background is observed which 
represents the total amount of the reaction pp AK 
including (j)(j) and (j)K^K~ contributions. These distri- 
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FIG. 9. Variation of the normalization integrals for the 
amplitudes describing the (^K'^ K~ (a) and (\><^ (b) amplitudes 
as functions of the incident p momentum. 
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FIG. 10. Combinatorial m(KK) distribution (six entries 
per event) for incident p momentum (a) below and (b) above 
1.5 GeV/c. The white region represents the (/>0 contribution, 
the grey region shows the (j)K^ K~ contribution and the black 
region shows the phase space (background -I- AK) contribu- 
tion. 



butions have been fitted using a Monte-Carlo generated 
shape for the iS.E distribution which includes simulations 
at all energies of ((xf) and AK final states, and a back- 
ground parametrization using the sum of two Gaussians. 
The number of non-resonant AK events was therefore ob- 
tained by 



4K 



N. 



<t>KK 



In the above expression, the number of AK events is the 
difference between the total number of events in each bin 
(Nt), the (N^^) and (j)K+K~ {N^kk) yields which 
were obtained from the channel likelihood fits, and the 
background {Nb) is drawn from the fits to the AS distri- 
butions. Due to the uncertainty in the background sub- 
traction, a 50% systematic error has been added quadrat- 
ically to the statistical errors. The background below the 
AK signal is relatively small, being in average of the order 
of 10 % increasing to 20 % only in the higher momentum 
regions. 
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FIG. 12. Cross sections in /ib for the reactions (a) 
pp — > 0<^, (b) pp — > K~ and (c) pp ^ 2K^2K'' cor- 
rected for unseen decay modes. 

VII. CROSS SECTIONS 

Having determined the number of events for each chan- 
nel we have computed the corresponding cross sections 
as: 

Events 



Acceptance x Luminosity 

Due to decreasing performance of the threshold Cerenkov 
counters in the last period of the data taking, part of data 
suffer an additional 20% systematic normalization error. 
These data represent about 30% of the total and have 
not been used in the calculation of the cross sections. 
However, properly scaled, these data can be used in the 
study of the angular distributions. These cross sections 
have been corrected for unseen (j) decay modes and are 
displayed in Table 1 and shown in Figs. [l^ . 

The (j)(f> cross section has been corrected assuming 
phase space in the calculation of the acceptance. This 
is not really a strong assumption as it can be seen from 
Fig. ^(c). In addition, a spin parity analysis of the (fxj) 
final state has been performed |^^. This analysis shows 
that the 00 system is dominated by J^'-^ = 2++. Cor- 
recting the mass spectrum with the results from the spin- 



TABLE II. (j>(j) cross 


section for the fine 


scan. 


cm. energy 


a(4>(f>) 


cm. energy 




GeV/c^ 


fib 


OeV /c 


fib 


2.219 


4.15 ± 0.34 


2.231 


3.73 ± 0.35 


2.221 


4.16 ± 0.36 


2.231 


2.97 ± 0.33 


2.221 


3.87 ± 0.54 


2.233 


4.20 ± 0.35 


2.222 


3.94 ± 0.35 


2.235 


3.52 ± 0.34 


2.224 


4.40 ± 0.36 


2.236 


2.48 ± 0.30 


2.226 


3.90 ± 0.34 


2.242 


3.38 ± 0.45 


2.226 


3.68 ± 0.34 


2.247 


2.68 ± 0.39 


2.228 


3.00 ± 0.30 


2.254 


2.63 ± 0.35 


2.229 


3.20 ± 0.32 







parity analysis has little influence on the shape of the 
integrated acceptance as a function of the 00 mass. 
Notice that: 

• The 00, (j)K^K^ and iK^ cross sections have dif- 
ferent shapes. The 00 cross section, in particu- 
lar, has a strong threshold enhancement while the 
(pK'^K^ and AK^ cross sections have a smooth in- 
crease as a function of the centre of mass energy. 

• The 00 cross section is rather large, about 3.5/i6 in 
the threshold region. 

• No evidence for narrow structures is found. 

The large 00 production close to threshold can be in- 
terpreted as a violation of the OZI rule. If the OZI rule 
is interpreted to forbid strangeness production in pp an- 
nihilations then the process pp —^ can proceed only 
via the small uu, dd component present in the wave 
function. With a deviation from ideal mixing 9 ~ 9q oi 
only 1° to 4", the is nearly 100 % ss. We can there- 
fore derive an upper limit tan^(6' — 9q) k, 2.5 x 10"^ for 
the ratio of cross sections a,^,^la^^ for production in pp 
annihilation. Although the cross section of pp louo has 
not been measured directly, an estimate can be obtained 



from the total 27r+27r 27r^ cross section 



which was 



measured to be about 5 mb in the energy range of our 
experiment. There are many reaction channels that con- 
tribute to this final state. If we estimate that it is 10 
% LJOJ then the expected 00 cross section is 10 nb, two 
orders of magnitude lower than our measurements. 



VIII. SEARCH FOR 5//,/ (2230) 

To search for the ^//,7(2230) resonance, a fine scan of 
the 00 cross section was performed during two different 
data taking periods. The results from these scans are 
shown in Fig. |l^ and displayed in Table II. 

Here the empty and full dots distinguish the two differ- 
ent sets of the data, showing good agreement in the size 
of the 00 cross section. No narrow structure is visible 
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TABLE I. Cross sections for pp (jxj), (j)K* K and 47^ 



Momentum region 


cm. 


Luminosity 




a((t)K+K') 


cm. 






Energy 








Energy 




GeV/c 


GeV/c^ 


nb"^ 




fih 


GeV/c^ 


/lb 


1.188-1.200 


2.147 


26.8 


2.86 ± 0.46 


0.25 ± 0.49 


2.158 


0.15 ± 0.15 


1.237-1.278 


2.168 


32.8 


3.45 ± 0.31 


0.96 ± 0.32 






1.300-1.330 


2.190 


69.1 


3.54 ± 0.18 


0.23 ± 0.21 


2.195 


0.56 ± 0.13 


1.360 


2.205 


33.4 


3.84 ± 0.24 


0.27 ± 0.26 






1.390-1.400 


2.218 


60.7 


3.55 ± 0.17 


0.63 ± 0.22 


2.219 


0.50 ± 0.17 


1.404-1.405 


2.220 


53.9 


3.88 ± 0.17 


0.54 ± 0.21 






1.410-1.415 


2.223 


38.0 


4.00 ± 0.20 


0.39 ± 0.25 


2.225 


0.52 ± 0.15 


1.420 


2.226 


25.7 


3.84 ± 0.24 


0.71 ± 0.29 






1.425-1.430 


2.228 


26.4 


3.09 ± 0.22 


0.97 ± 0.32 






1.435 


2.231 


47.7 


3.33 ± 0.17 


1.16 ± 0.235 


2.233 


0.58 ± 0.15 


1.440-1.450 


2.235 


35.1 


3.41 ± 0.19 


1.19 ± 0.26 






1.465-1.480 


2.242 


48.7 


3.02 ± 0.15 


1.07 ± 0.22 


2.247 


0.44 ± 0.26 


1.500-1.506 


2.255 


42.7 


3.01 ± 0.16 


1.77 ± 0.25 






1.550 


2.272 


24.3 


2.25 ± 0.17 


1.90 ± 0.31 


2.282 


0.34 ± 0.33 


1.600 


2.289 


20.1 


2.15 ± 0.21 


2.40 ± 0.40 






1.650 


2.307 


13.6 


1.97 ± 0.24 


3.26 ± 0.45 






1.700-1.750 


2.333 


32.6 


2.04 ± 0.14 


3.56 ± 0.30 


2.348 


0.95 ± 0.52 


1.800 


2.360 


39.3 


2.04 ± 0.12 


4.30 ± 0.25 






1.900-1.950 


2.404 


15.4 


1.61 ± 0.18 


3.63 ± 0.41 


2.424 


2.55 ± 0.93 


2.000 


2.430 


52.1 


1.93 ± 0.11 


3.91 ± 0.24 








2.26 



mC* <P) GeV/c' 



in the data. We have fitted the (fxj) mass spectrum using 
a polynomial and a Breit-Wigner form representing the 
^//j(2230) with m=2235 MeV and r=15 MeV, param- 
eters measured in the reaction J/4' —>■ ^S. where £, ^ pp 
by the BES experiment [^. The Breit-Wigner form is 
written as: 



<7bW = (WiWf) X 



(2J+1) 



(25i- 



1)(2S2+1) 

p2 



Q- — zLm"^ 



Here (wiWf) is the double branching ratio, WiWf — 
BR(X pp) X BR(X (jxj)). The Si terms are the spins 
of the initial proton and antiproton (1/2), and J is the to- 
tal angular momentum of the resonance, reducing the an- 
gular momentum term, (2J-|- 1)/ ((2S'i + 1)(2S'2 + 1)) to 
5/4 in the case of a J = 2 resonance. A limit on the prod- 
uct of the branching ratios of: BR(^ pp) x BR(^ — > 
M) <&x 10-5 at 95 % c.l. is obtained. 



IX. CONCLUSIONS 



FIG. 13. Cross section in /ib for the reaction pp (fxj) 
in a fine scan over two different periods of data collection 
corrected for unseen <j) decay modes. Open circles: 1991 data; 
full circles: 1993 data. The line is the result from the fit 
described in the text, the curve represents a Breit-Wigner 
resonance whose amplitude is at the 95 % c.l. upper limit for 
the production of a ^//j(2230) with a mass of 2235 MeV and 
a width of 15 MeV. 



We have performed a high statistics study of the reac- 
tion pp AK^ using in-flight p from 1.1 to 2.0 GeV/c 
incident momentum interacting on a hydrogen jet target 
of the JETSET (PS202) experiment at CERN/LEAR. 
The reaction is dominated by a strong 00 production at 
threshold whose strength exceeds by two orders of mag- 
nitude the yield extracted from a simple interpretation 
of the OZI rule. 
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Several models have been proposed in order to explain 
the large OZI violations observed in hadron induced re- 
actions and particularly in some pp annihilation chan- 
nels. However, few quantitative calculations exist for the 
specific channel under investigation in the present exper- 
iment. 

A model which interprets pp annihilations to (ficj) as 
due to KK rescattering |^ is able to predict the or- 
der of magnitude of the cross section (« 2.4/i6), but 
not the detailed shape of the observed spectrum. Other 
models make use of hyperon-antihyperon intermediate 
states In this case the size of the cross section 

is underestimated by a factor of about 4. Other ways 
to enhance production of the (fxj) system have been sug- 
gested in ref. jl^] by invoking the hypothesis of intrinsic 
strangeness content in the proton. The authors suggest 
that (/)(/) production could originate from rearrangement 
diagrams with strange quarks originating from the pro- 
ton sea which are polarized with total spin S=l. The 
authors conclude that these connected rearrangement di- 
agrams very likely mask any possible glueball resonance 
contributions which are expected to be dominant among 
the disconnected diagrams. Further information and pos- 
sible new inputs to the problem may come from a spin 
analysis of the observed 00 threshold enhancement . 

No evidence for narrow resonance is found and we set 
an upper limit for the production of ^//,7(2230) with 
m=2235 MeV and r=15 MeV of: BR(^ ^ pp) x BR(^ 

< 6 X 10~^ at 95 % c.l. Combined with our scan of 
pp KgKg a consistent and stringent rejection of 
a narrow resonance (F < 30 MeV) appears. 
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